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Abstract 

The sine-Gordon equation in light cone coordinates is solved when 
Dirichlet conditions on the L-shape boundaries of the strip {t G [0, T]} U 
{x G [0,oo]} are prescribed in a class of functions that vanish (mod 2-k) 
as x — > 00 at initial time. The method is based on the inverse spectral 
transform (1ST) for the Schrodinger spectral problem on the semi-line 
x > solved as a Hilbert boundary value problem. Contrarily to what 
occurs when using the Zakharov-Shabat eigenvalue problem, the spectral 
transform is regular and in particular the discrete spectrum contains a 
finite number of eigenvalues (and no accumulation point). 



1 Introduction 

The sine-Gordon equation in light cone coordinates (dimensionless variables, 
subscript stands for partial derivative) for the field 9(x, t), 

xt + sin0 = O, (1.1) 

is one of the most widely studied nonlinear wave equation because of its intrinsic 
mathematical properties and its wide applicability in physics. The initial value 
problem on the infinite line x G (— 00, +00), namely the initial datum 9(x,0) 
together with the asymptotic behaviors 9 — 2mr as x — > ±00, has been solved 
by the inverse spectral transform (1ST) in jj], 0]. 

The boundary value problem for ( |l . 1| ) on the quadrant {x > 0, t > 0} 
has been solved in Q and reconsidered in In both approaches it has been 
shown that a constant boundary datum 9(0, t) on an initial vacuum 9(x, 0) = 
produces an infinite series of discrete eigenvalues with accumulation point at 
the origin. The problem is then that the Jost solutions are not continuous in a 
neighborhood of the real axis and the spectral theory fails. 
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We solve this problem by considering an alternative Lax pair for the sine- 
Gordon equation deduced from the work || on two-dimensional extensions of 
sine-Gordon. The principal spectral operator is the Schrodinger eigenvalue prob- 
lem with complex potential that we reconsider on the half-line x > with a given 
boundary value of the potential in x = 0. 

The spectral theory developped here follows the approach introduced in || 
for stimulated Raman scattering, reformulated in j| , and later applied to others 
continuous @ or discrete |§| coupled wave systems. In the next section the Lax 
pair is set and the boundary value problem is discussed. Section || is devoted to 
the solution of the direct spectral problem, namely to the method of construction 
of the spectral transform, or nonlinear Fourier transform, of the sine-Gordon 
field. 

The time dependence of the spectral transform is then derived in section || 
where it appears to be a nonlinear Riccati type evolution. It is examined for 
particular boundary data and shown to lead asymptotically to a finite set of 
discrete eigenvalues (soliton component of the spectrum) . In section ^| we solve 
the inverse problem which is the reconstruction of the sine-Gordon field from 
the spectral transform. This is done by formulating the reconstruction as a 
Hilbert boundary value problem. In doing this attention must be paid to the 
asymptotic boundary values in the spectral plane required for the solutions of 
this Hilbert boundary value problem. Indeed one of the Jost solution has an 
essential singularity for large real values of the spectral parameter. Still the 
solution can be obtained by a Cauchy-Green integral equation with a contour of 
integration that differs from the one used in the infinite line case. The spectral 
theory of the Schrodinger operator on the half- line can also be developped within 
the GePFand-Levitan-Marchenko fromalism as shown in Q where the boundary 
value problem for the Korteveg-de Vries equation is studied. 



2 Alternative Lax pair for sine-Gordon 

Weakly commuting operators. Let us consider the following pair of scalar 
operators: the principal spectral operator 

T 1= d xx + k 2 -p, (2.1) 
where p(x, t) is the principal potential, and the auxiliary spectral operator 

T 2 = d xt + ad t + P , (2.2) 
where a(x, t) and [3(x, t) are the auxiliary potentials. 

Originaly introduced for two-dimensional problems |lC| |, the weak compati- 
bility of those operators results in requiring that they commute on their common 
eigenspace. In other words 

[T 1 ,T 2 ] = fT 1 +gT 2 , (2.3) 
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where / and g are two unknown functions. Straightforward examination of 
highest order derivatives gives 

/ = , g = 2a x , (2.4) 

hence we are left with the operatorial equation 

[T 1 ,T 2 ] = 2a x T 2 . (2.5) 

This equation exlpicitely implies the following system 

Px = -a xx + 2aa x , (2.6) 
Pt = -2Px , (2-7) 
Pxt + apt = -/3 XX + 2a x (3 , (2.8) 

in which two equations allow to eliminate two fields and the remaining equation 
is the nonlinear evolution. 

The sine-Gordon equation. Assuming for all t the asymptotic behaviors 
(we will prove later that these behaviors result from the properties of the initial 
datum) 

oo : {p^O.a-tO.aj-fOJ^O}, (2.9) 
two of the above equation can be integrated and give 

/9 B = -2a/3, p = a 2 -a x , (2.10) 
while the remaining equation reads as the evolution 

Pt - ~2(3 X . (2.11) 
To eliminate two fields we define the new variable 6(x,t) by 

P=\e w , (2.12) 



which from (2.9) will be assumed to obey 

6 X — » . (2.13) 



The relation ( p. 10 ) then gives p and a in terms of 9 as 

i „ i „ 1 



p2 



(2.14) 



where one can check that the behavior (2.13) ensures the asymptotic expressions 



The evolution p t — —2f3 x can then be factorized as 



{ ~L + ie *}( 6 *t + sin6 ) =° • ( 2 - 15 ) 

This can be integrated once and the integration constant can be evaluated by 
reading the evolution as x — > oo where 9 — > A. The results is the following 
complex extension of the sine-Gordon equation 

9 xt + e iX sin(0 - A) = , A = lim 9 (2.16) 



associated to the asymptotic behavior ( 2.13 ). Note that in general the parameter 
A is a function of time t. 

In the case of real valued 9, and hence real valued A, the above evolution 
implies A = and hence reduces to 

9 e K => 6 xt + sm6 = , 9 — ► mod 2tt . (2.17) 

x — >oo 

Thus the standard sine-Gordon equation results as the weak compatibility of 
operators T\ and T2 for real valued field 9. 

It is important to remark here that the imposed vanishing asymptotic be- 
havior of the derivative 9 X for all t as given in ( [2.13 ) actually results from the 



above assumed asymptotic behavior of 9 as soon as 9 X vanishes at infinity at 
initial time. Indeed, by integration of the sine-Gordon equation we can write 



9 x (x,t) =6 x (x,0) - [ dr sin 9{x,t) . (2.18) 
Jo 

This expression serves actually as a means to go from Dirichlet to Neuman 
conditions on the boundaries x — and x = 00. It is also useful to prove that 
the assumed vanishing asymptotic behavior of the principal potential p(x, t) 



results from the asymptotic value of 9{x,t). Indeed, upon derivation of (2.18) 
it is easy to obtain that if the initial datum is such that 9 x (x, 0) and 9 xx (x, 0) 
vanish as x — > 00, then 9 XX vanishes for all time as x — > 00. 



Note also that the expression (2.16) of the sine-Gordon equation for complex 



variable can be put in the form (2.17) by the following change of variables 

rt 



9- A , t = / dt'e^ 1 '^ , (2.19) 







which gives 



(p XT + simp = , <p — > mod 2tt . (2.20) 



We have assumed hereafter the general situation where the constant A is allowed 
to depend on time t. 
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Dirichlet boundary data. We consider then the sine-Gordon equation (1.1), 
a priori for a complex field 9, on the strip 



{xe [0,oo]}u{ie [0,T]} , 



(2.21) 



for any value of T. The Dirichlet conditions are the data of the function 9(x, t) 
on the L-shape boundaries of this strip, namely 



t = Q : 6(x,0)=6(x), x = : 0(0, t) = d (t) 



(2.22) 



The auxiliary spectral operator ( |2.2| ) explicitely contains the variable a(cc, i) 
which has then to be evalu ated in x — explicitely in terms of the Dirichlet 
data ( 2.22| ). By expressing ( [2.18 ) in x = we have 



:(0,t) =0 X (O) 



dr sin 9q (t) 



(2.23) 



which from (2.14) gives 



«o(*) = = / dr sinfloM - -9 X (0) 



(2.24) 



where ao = a(0, t) is then completely determined from (2.22). Similarly we have 
to evaluate po(t) = p(0, t) occuring in the principal Lax operator. By derivation 
of (2.18) we obtain 



Po(t) = ^xx(O) - - / dr 9 x (0,t)cos9 (t) - ~^(0,*) , 



(2.25) 



where 9 X (0, t) is explicitely given in ( 2.25 ). We have then obtained that the 
data of 9q (t) and 9(x) determine completely the boundary values in x = of 
the potentials of the fundamental operators T\ and T^. 

For the other boundary x = oo we first select a class of functions for the 
initial datum 9(x). As usual for implementation of the spectral transform, let 
us assume that 9(x) belongs to the Schwartz space S (bounded functions that 
vanish with all their derivatives faster than any power of 1/x). Although this 
condition is rather stringent, it fits well a generic physical situation, in particular 
the case of an initial vacuum. It is then easy to obtain from ( [2.18] ) the property 



{9(x)eS, 0(oo,*) = O mod 2tt| => 6^(00,^=0 



(2.26) 



This can be worked again for the second derivative and finally the potentials 
p(x,t) and a(x,t) do vanish as x — > 00 and (3 goes asymptotically to 1/4. 
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3 Spectral problem 



The Schrodinger spectral problem on the half- line Titp = with complex po- 
tential, namely 

x>0, ip xx + {k 2 - p)*P = , p(0)=po, (3.1) 

is solved now and the inverse problem will be takled by means of a Hilbert 
boundary value problem. A partial solution has already been given in in the 
case of real valued potential. We rework everything here with slightly different 
definitions and in particular give the complete solution of the Hilbert boundary 
value problem. Note that this spectral problem has also been considered in . 



Basic solutions. We redefine the spectral equation as 

i> = e lkx , <f> xx + 2ikcj) x - p<j) = , (3.2) 
and select two basic solutions </> + and <f>~ by 

^-(k 1 x) = l + ±-J*dy<j ) -(k 1 y)p(y)[l-e- 2 ^ x -y^ , (3.3) 



dy (/) + (k 1 y)p(y) 



1 f°° 

+ 2ikJ dy <t> + (^y)p{yy~ 2ik{x ~ v) ■ ( 3 - 4 ) 

where conditions on p{x) are assumed (e.g. p(x) belongs to the Schwartz space) 
that ensure uniqueness of the solution of the above integral equations. It is 
useful to compute the derivatives 



4>-{k,x)= dy^-(k, y )p(y)e- 2ik(x -y^, (3.5) 
Jo 

POO 

<j>+(h,x) = - dy cb+ (k, y)p{y)e- 2lk ^ . (3.6) 



An important property for the following is that another set of basic solutions 
is constituted by ^(—k, x)e~ 2lkx . Another important property is the fact that 
the wronskian of two solutions of (3.1) ia a constant, which means for any two 
solutions (f>i and 2 of (3.2) 

= 0m02 - <t>ifo,x = C e- 2lkx , (3.7) 

where C is a constant. 



Spectral transform. The solution <fi + is called the physical solution and 
serves to define the scattering coefficients as follows. The boundary value in 
x = 

x = : = l + p(fc) , (3.8) 



6 



defines the reflection coefficient p(k) as 

1 



p(k) 



2ik 



dy P (y)<j> + (k,y)e 2ik y 



and consequently (3.6) gives 

x = : 0+ = -2ikp{k) . 
The boundary value in x = oo 

x — > oo : + — > r(fc) , 
defines the transmission coefficient r(k) as 



1 f 00 



and we have also 



x — ► oo : 0^ — > 



(3.9) 

(3.10) 
(3.11) 

(3.12) 

(3.13) 

2ikx 



The wronskian relation (3.7) applied to the solutions </> + (/c) and <^> + (— fc)e 
at x = and a; = oo successively gives 

p{k)p(-k) + r(fc)r(-fc) = 1 , (3.14) 

called the unitarity relation. 

From its very definition, the spectral tranfrom p{k) seems to be ill defined 
in k = 0. Actually we recall in the appendix [A| that 



k = : p= -1 , r = 0, 
which is well known in the infinite line case fll2l . 



(3.15) 



Hilbert boundary value problem The solution $T is called the intermedi- 
ate solution and serves to complete the basis. We can for instance express <j) + {k) 
in terms of the basis 4>~(k) and 4>~(— k)e~ 2tkx , the coefficients being computed 
by comparing the values in x — of both the functions and their derivatives. 
We obtain 

(/} + (k ) x)-(f>-(k,x) = p(k)e- 2lkx <j)-(-k,x) , (3.16) 

which constitutes a Hilbert problem of the real fc-axis, we shall come to that 
argument later. For the moment we express the above relation by changing k 
to — k and solve the resulting linear system for the function 4>~ {k) to obtain 

cf>-(k,x) = _ 1 [4>+(k,x)-p(k)e- 2ikx 4>+(-k,x)] . (3.17) 
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Using the definitions of <f> and <f> x and the above relation we can compute the 
different boundary values and get 



x = 



1 



-(-k) r(k) 



= , 



T{k) 



(3.18) 
(3.19) 



where we have also used the unitarity relation (3.14). 



4 Solution of sine-Gordon 

We have constructed in the preceding section the spectral transform p(k) of the 
principal potential p(x). We shall prove in section ^| that the datum of p(k) in 
the upper half plane completely determines the potential p(x), and we shall give 
the tools to reconstruct this potential. But we consider first the question of the 
time dependence of the spectral transform when the potential depends on time 
according to the evolutions ( 2.1CQ , reducing to sine-Gordon in the real valued 
variable 9. 



Evolutio n of the spectral transform. From the weak compatibility com- 
mutation (2.5) of operators T\ and T 2 , the selected eigenfunctions f/' ± are not 
eigenfunctions of T 2 and our purpose here is to determine their proper time 
evolution. Equation (2.5) implies the property 



Tyip — => (Ti - 2a x )T 2 ip = , 
and thus using the remarkable identity 

T x i(i = => (Ti - 2a x )(ip x + a^j) = , 



(4.1) 



(4.2) 



which simply results from the structure of T\ in terms of the potential a, we 
may write 



T 2 ^ + = n (V>++cn/; + ) 



(4.3) 



Hereabove, the ^-independent coefficient Q(k, t) is called the generalized dis- 
persion relation. Note that one should actually expand T2ip + on a basis, say 
ipx(k) + aip + (k) and ip+(—k) + aijj + {— k), but examination of the asymptotic 
behaviors as x — > oo readily reduces the expansion to the form ( [1.3] ) and gives 
the evolution 



Ti 



i 

Aik 



(4.4) 



of the transmission coefficient. 
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The value of the dispersion relation fl together with the time evolution of 
p are obtained by expressing equation ( |4.3| ) and its x-derivative in x = 0. We 
obtain 

(a ~ ik)p t +Po(l + p) = n[oo(l + p) + ik(l ~ p)] , (4.5) 
(og - ika - k 2 )p t +2a (3 (l + p) + ik/3 Q (l - p) = 

n[(a 2 -k 2 )(l+p)+ika (l-p)} . (4.6) 

where we have used the expressions p — a 2 — a x and j3 x — —2a/3 and the 
equation Tiip = 0. This is a system of equations for the unknowns {pt, 51} in 
terms of the boundary data ao and @q and of p. Note that it is compatible with 
(see appendix [A|) 

k = {p = -l, t = 0} . (4.7) 
Multiplying the first equation by ao and substracting with the second gives 



a 



pt = (l+p) V + Po^ -tt(I-p) • (4-8) 



ik 



Replaced in any of the preceding equations, we obtain the dispersion relation 
which after some algebra can be written 

n = p ^ a °~ 2lk) + it> l+2k2) - (4 - 9) 

Thus the evolution of the reflection coefficient eventually reads as the Riccati 
differential equation 

p t = ap 2 +bp + c , (4.10) 
with entries which depend on the boundary data by 

a = 2iP" ( ° ~ ) ' 6= ^p (a o + 2A:), c = -^p- (a + 2tfe) . (4.11) 



Soliton generation. To illustrate the above formalism we consider the case 
of a piecewise constant boundary datum, i.e. 



t = 



(x) = o , t < h 



)(t) = , t > tx ■ e (t) = v 



(4.12) 



For such a simple case, the Zakharov-Shabat spectral transform has been shown 
to generate an infinite number of discrete eigenvalues with k — as an accumu- 
lation point jy . 



Here we have from (2.24) 



t < h : a (t) = 0, /3 = - , 



t > ti : ao(t) = — (i — ti) sin 99, /3q 



(4.13) 
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For such expressions, the time-evolution ( [1.10 ), with vanishing initial condition 



p(k, 0) = 0, cannot be explicitely solved. However we can evaluate the large 
time asymptotic behavior as indeed 

b _ al + 2fc 2 



= 2- 

xq(oco — 'lik) 

^ -4- 9ik 

1 . (4.14) 



a ao(ao — 2ik) t 
c olq + 2ik 



a o;o — 2ik t^oc 

The resulting asymptotic evolution then reads 

p t ~a(p+lf, (4.15) 

which is easily solved to eventually get for large time t 

^ -e^siny[|(t-ti)siny-fc] {t - hf 
P[ ' ' ~ 32ik 3 + e i <Psmtp[l(t-t 1 ) sin <p-k] (t - <i) 2 ' 

This function has 3 simple poles given by 

= 32iZ^ - 2e^sin<p Z„ + ie^sin 2 ^ = . (4.17) 

o 

For different values of ip, one or two of these poles are in the upper half plane 
corresponding to a solution with a one-soliton or a two-soliton component. 



Note that the Riccati time evolution ( 4.10 ) is not in general explicitely solv- 
able. Indeed, simple explicit solution could only be obtained for b/a and c/a 
time independent which would immediately lead to 8 (t) = 0. Such does not oc- 
cur in the Zakharov-Shabat spectral transform of sine-Gordon where piecewise 
constant boundary data in x = produce an explicitely solvable evolution Q. 
Then it has been shown that the spectral transform contains an infinite series 
of poles, accordingly with see also [[l3). The problem is that these poles 
possess an accumulation point (the origin k = 0) and hence some of the basic 
solutions are not continuous in the vicinity of the origin. Using the present spec- 
tral transform based on the Schrodinger spectral problem removes this problem 
and thus provides the correct solution to the Dirichlet boundary value problem 
for sine-Gordon in light-cone coordinates. 



5 Inverse problem 

Solving the inverse problem will consist in computing the potential p(x) from the 
datum of p{k). This is done by solving first the Hubert boundary value problem 



(3.16) which actually requires the knowledge of the analytical properties of the 
functions ^ ± (fc) together with their behaviors at large k. 
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Analytical properties. The Volterra integral equation ( |3.3| ) allows to deduce 
that 4>~ is an entire function of k bounded in the lower half plane. The function 
<f> + is not analytic but the function cf) + — (f> + /t solves the Volterra integral 
equation 

0+ = 1 - — J dy <t>+{k, y)p(y) [l - e - 2 «(-rtj . (5.1) 

Hence <j) + /r is analytic in the upper half-plane and thus, from definition ( |3.12| ). 
1/t is also analytic. Consequently cf> + has a pole where 1/r has a zero, these 
are the discrete eigenvalues of the spectral transform named k n . 



The definition (3i)) of the reflection coefficient shows that p can be extended 
to the upper half plane (unlike in the infinite line case where the integral would 
run on the whole real axis) where it possess the discrete eigenvalues as poles. 
The integral equation (3.4) can be rewritten as 

</>+ - 1 + pe- 2lkx + f dy cj>+(k, y)p{y)[l - e -**C*-v>] . (5.2) 



2ik _ 

By taking the residue of the above in k n we get 

e 2 ^Res(0+) = Pn - -L f X [e 2 ^yRes(cb+)}p(y)[l - e 2ik ^~^] , (5.3) 



where obviously 



Res(p) . (5.4) 



By comparaison of the above integral equation with the one for the function 
4>~(— k,x) read in k = k n we immediately conclude that 

e^Res^+ik.x)) = p n <f>- (-k n ,x) . (5.5) 

fen 

This closes the analytical properties needed to the solution of the Hilbert bound- 
ary value problem. 



Asymptotic behaviors. It is necessary now to determine the large k behav- 
iors of the solutions ± . By repeated integration by part, the solutions obey 



^~ 1 + 27^ (1) + W 



(2) 



(2) _|_ (p (2) e -2 l fc a: 



+ 



with the following coefficients 



s(2) 



dy p(y) , 

p(x) + dy p(y) I dz p(z) 



<^ (2) = Pa 



(5.6) 
(5.7) 

(5.8) 

(5.9) 
(5.10) 
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The solution of the Hilbert problem will have to verify those behaviors, with in 
particular the essential singularity for (j)~ on the real axis. 



The relation (5.S) furnishes the means to obtain the potential p from the 
solution <j) + namely 

p(x) =d x 4> {1) {x) , (5.11) 

where is the coefficient of l/2ik in the large k expansion of cf> + (k,x). Note 
finally that the spectral transform p(k) obeys in the upper half-plane 

m = -Wr + ° {l/k5) ' (5 ' 12) 



obtained again by integration by part on the definition (3.9) with help of the 
behavior fl5.6| ). Note that this is consistent with the boundary value of <j) + in 
x = given in (ED3) and the large k behavior fl5.6|). Indeed 



>l+ 1 _ Po , 

1 + \<t> {l) + — L^ (a) + >l P ° 



fc->oo 2ik (2ik) 2 x^o (2ik) 2 

This works as well for 0~ , hence in summary, the limits k — > oo and 
commute. 



Solution of t he Hi lbert problem. The solution (jr 1 of the Hilbert boundary 
value problem (|3. 16 ) which verifies the preceding analytical properties and large 



k behaviors is given from the data of p(k) in Im(fc) > (including eigenvalues k n 
in the upper half-plane and related residues p n ) by the Cauchy-Green integral 
equation 

(f) -(k,x) = l + — frf (-A, x)p(\)e- 2 * x * 

2nr J-oo+io X - k 



N 



-2ik n x 



+ J2 , , 4>-{-K,x), lm(fc)<0, (5.13) 

n—1 

and the explicit formula 

-i p+oo-M ii 

^{ k ,x) = l + — <t)~ (~\, x)p(\)e~ 2iXx 

2*71" J-oo-tO X - k 



Y, , —r <t>-(-k n ,x) , lm(fc)>0,. (5.14) 

^ ^ tin & 



As shown for the first time in 0, the choice of the particular contour in the 
complex plane for the above expressions is essential. 
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Indeed, while the factor e~ 2lXx produces a bounded kernel for x > in the 
lower half A-plane, it does not so in the upper half A-plane, that is for the 
expression of <j)~ . Assuming no pole for p on the real axis, we can use contour 



integration on (5.13) to write 



(f>~(k,x) = 1 + 



-2iki 



dX 



<t> (-X,x)p(X)e 



-2i\x 



<f> {—k,x)p{k)e 



E 



k n k 



cf> (-k n ,x). 



(5.15) 



The first consequence is that indeed these solution obey the Hilbert relation 
( |3.16l). T he second consequence is that they do verify the asymptotic behav- 



iors (5J3) and (5.7) obt ained here simply by expanding the ratio 1 /(A — k) in 
expressions ( 5.14 ) and ( 5.15 ) where kernels are bounded, unlike in ( 5.13j ). 



The analytical properties naturaly follow the expressions ( |5.14|) and ( |5.13|) . 
Indeed, as Im(fc„) > 0, the analyticity of the solution <fi~(k,x) of ( 5.13 ) simply 
result from the property (continuity, boundcdness) of p(k) on t he real axis. 
Then the solution (j> + (k,x), given by the explicit expression (5.14), appears as 
the sum of an analytic function and the contribution of its simple poles k n . 
Finally the relati on (|5.5 [) between the residue of 4> + in k n and p n readily follows 



from expression (5.14) 



Reconstruction of 9(x,t). Once the spectral transform p(k,t) has been ob- 



tained by solving the Riccati time evolution (4.1C) for the initial datum p(fc,0), 
the solution of sine-Gordon 8(x,t) is built out of the integral equation (5.13) 
as follows. The relation ( |5.11 ) gives the value of p(x,t) from the coefficient 



(j^ x \x, t) of \/2ik in the large k expansion of 4> + {k, x, t), namely 

f+oc — iO 



</> (1) (x, t) = / d\ 4>~ (-A, x, t)p{X, t)e 



2i\i 

oo— iO 



N 

-2iY,^(~k n ,x,t) Pn e- 2ik " x . (5.16) 

n=l 

This expression shows the two building components of the solution: the radia- 
tion component represented by the continuous integral and the soliton compo- 
nent represented by the discrete sum. Note that the number N and positions k n 



of the solitons depend on time as allowed by the Riccati time evolution fl4.10|) . 
The simplest way to obtain now the solution 8(x, t) is to use the relations 

Pt = -2f3 x , p=4 1 >, (5.17) 
easily integrated with the values in x = of [3. We eventually arrive at 

e ifl(*,t) = e ifla(t) _ 2d t <t> {1) {x,t) . (5.18) 
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This provides then the solution 9(x,t) from the boundary value 8o(t) and the 
solution 4>~(k, x, t) of the integral equation (5.13) completely determined by 
p(k,t). 



6 Conclusion 

The solution of the Dirichlct boundary value problem for the sine-Gordon equa- 



tion has been shown to be completely characterized by the time-evolution (4.10) 
of the spectral transform p(k, t) evolving explicitely the boundary datum 8(0, t), 
with an initial datum p(k,0) determined from the initial condition 9(x,0). 

This solution has many important aspects. First the solution is stable for 
initial conditions in the Schwartz space, much like what occurs for the infinite 
line case. Second, the time evolution of the spectral transform is regular in the 
sense that it does not produces accumulation of poles. Last but not least this 
approach allows to treat the finite line case as indeed it is an open-end boundary 
value problem. This property was already remarked in Q and later discussed in 
[pi , and mainly results from the fact that the solution is completely determined 
by the integral equation (5.13) where the x-dependence is parametric. 

In order now to descibe in more details interesting generic behaviors of the 
solution of sine-Gordon, we would need to study the Riccati time evolution 



(4.10) where the coefficients are time-dependent. This aspect of the work is 



reported to further studies as it would require use of numerical simulations. 



A Behaviors at the origin 



n( x ) = 1 + Pa - / dy(x - y)4>t (y)p(y) 

J X 

4>t( x ) = [ 1 + po\4>o ( x ) > 



The definition (3.4) of the solution <p + (k,x) and the Hilbert relation (3.16) 
written both in x = lead to the system (in k = we note 4>(0, x) = ^o(^) and 
P(0) = P0 ) 



(A.1) 
(A.2) 

(A.3) 
does not 

(A.4) 



while (f> solves 



(f> (x) = 1 + / dy(x - y)(/) (y)p(y) . 
Jo 

As the above Volterra integral equation has a unique solution, 
vanish and the equation (A.l) can be written by using (A.2) as 

/>OC 

0+(a;) [1 - l/0o (x)] = - / dy{x - y)cp+(y)p(y) . 
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This has the unique vanishing solution and consequently we have proved that 

4(aO=0, Po = -l- (A.5) 
Finally the unitarity relation (|3.14j ) leads immediately r(0) = as expected. 

B Reduction 



We determine here the symmetry properties of the basic solutions and the spec- 
tral transform resulting from the constraint 9 G K. First we note that 



6 G K T x * = Ti - 2^ , 
where for any function f(k) of the complex variable k, we define 

r = 7(k) . 



(B.l) 



(B.2) 



Second the operator T\ has the property ( |4.2| ) and consequently when 6* is real 
valued (or equivalently when a is purely imaginary), the eigcnfunctions ip* and 
•ipx + atp are linearly dependent, namely 



(d x + a)ip+(k) = Aij>{k) + BV (-*) 
(ft, + a)V~(fe) = ^(fe) + Djj + (-k) 



(B.3) 
(B.4) 



The coefficients hereabove are computed by expressing these expansions in x 
oo with help of (3.11) and ( 3.1 6| ) . We obtain using a(oc) = 



A = , 



C = ik- 



B = ik 



p(k) 



T(k)T(k) 



D = 



r(k)_ 
f(-fe) 
ik 



r(-fc)r(-fc) 



(B.5) 
(B.6) 



The next step consists in expressing (B.3) and (B.4) in x — which results 
in a set of two equations linking the spectral coefficients and the boundary 
value ciq. By eliminating «q between these two equations and after some simple 



algebra using the unitarity relation (3.14) we eventually obtain first 
l[l + p(k)][l+p(k)} , l[l + p(-fe)][l + /9 (-fc)] 



-(fc)r(fc) 



-(-fc)T(-fc) 



1 



(B.7) 



which constitutes the condition for a spectral transform p{k) to ensure a reaZ 
valued reconstructed field 9(x). The second equation gives a relation between 
the boundary value ao and the spectral transform, namely 



ao = —ik 



1 - 



T(-k)T{-k) 



p{k) 



1 +?(*-) " 
r(fc)T(fc). 



(B.8 
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This relation must in particular be valid as k — ► 00 where the behaviors of 
the spectral coefficients p and r are given by 



— H h O(-) 

2ik 2(2ik) 2 ( k 3) 



(B.9) 
(B.10) 



The behavior of r is obtained from its definition ( 3.12 ) an d the large /c-behavior 
( JTbI ) of the eigenfunction <j> + . By using then the relation (5.8) it is not difficult 
to obtain 



7 = «o + 



dx a 



(B.ll) 



The reduction constraint 9 £ R is here equivalent to 7 — 7 = 2ao • With all this 
in hands, the above relation (B.E) does hold as fe — + 00 where the zeroth order 
term in the right hand side is indeed oto. One can check then that the next 
order is a vanishing identity as we have within the reduction 



Pa + Pa - 2a a = 



(B.12) 
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